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Abstract 
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Introduction 

The possibility of constructing and (at least partially) solving by algebraic and/or analytical methods, 
one-dimensional interacting quantum spin chains, is one of the major achievements in the domain of 
quantum integrable systems. Its main tool is the quantum i?-matrix, obeying a cubic Yang-Baxter 
equation, the "coproduct" properties of which allow the building of an L-site transfer matrix with 
identical exchange relations and the subsequent derivation of quantum commuting Hamiltonians [1]. 
The same structure is instrumental in formulating the quantum inverse scattering procedure, initiated 
by the Leningrad school [2]. 

A subsequent development was the definition of exactly solvable open spin chains with non-trivial 
boundary conditions. These are characterised by a second object: the reflection matrix K, obeying 
a quadratic consistency equation with the R matrix, with the generic abstract form RKRK = 
KRKR [3-7]. Using again "coproduct-like" properties of this structure one constructs suitable 
transfer matrices yielding (local) commuting spin chain Hamiltonians by combining K and semi- 
tensor products of R [4]. 

Many efforts have been devoted to this issue [8-23], based on the pioneering approach of Sklyanin 
and we here aim at treating a particular, but very significant, class of examples for this problem. 

To better characterise the type of spin chain which we will be considering here it is important to 
recall that both R and K matrices have an interpretation in terms of diffusion theory for particle- 
like objects identified in several explicit cases with exact eigenstates of some quantum integrable 
field theories such as sine-Gordon [24], non-linear Schrodinger equation [4,25] or principal chiral 
model [26]. R describes the basic 2-body scattering amplitudes and K describes a 1-body scattering 
or reflection on a boundary. The Yang-Baxter equation (YBE) and reflection equation (RE) then 
characterise consistent factorisability of any L-body amplitude in terms of 1- and 2-body scattering 
amplitudes, regardless of the order of occurrence of the 1- and 2-body events in the diffusion process. 

As a consequence, when one describes the scattering theory of a model with more than one type 
of particle involved, one is led to introduce several operators of R and K type. The case which we 
examine here corresponds to a situation where the states involved can be split into resp. particles 
and anti-particles with a suitable representation of CP transformation acting on the states. Within 
the context of integrable field theories it is justified to denote them respectively "solitons" S and 
"antisolitons" A. Assuming that the 2-body diffusion conserves the soliton or antisoliton nature of 
the particles, but that the reflection may change it, one should therefore consider four types of R 
matrices (resp. -Rff, -Ra5; ^sai ^aa) connected by CP operations; and four reflection matrices 
(resp. i^f , Kj^, Kg, K^). It then becomes possible to deflne several non-equivalent constructions 
of commuting transfer matrices. As a consequence, one sees that a variety of spin chain models can 
be built using a Sklyanin-like procedure, depending on which transfer matrix is being constructed 
and which reflection matrices are used to build it. Locality arguments also come into play, leading 
to more complicated combinations of transfer matrices as we shall presently see. 

We shall here describe the construction, and present the resolution by analytical Bethe Ansatz 
methods [16,27], of open spin chains based on the simplest rational R matrix solutions of the Yang- 
Baxter equations for underlying sl{N') and sl{M.\N') Lie (super) algebras. These solutions, with a 



rational dependence on the spectral parameter, are instrumental in defining the Yangian [28]. We 
shall consider two types of associated reflection matrices K to build two distinct types of integrable 
spin chains: one which entails purely soliton- and antisoliton-preserving reflection amplitude by two 
matrices Kg and K^ (hereafter denoted "soliton-preserving case" or SP); the other which entails 
the two soliton-non-preserving reflection amplitudes Kg and K^^ (hereafter denoted "soliton-non- 
preserving case" or SNP). Closed spin chains based on sl{M\M) superalgebras were studied in 
e.g. [29] and, in the case of alternating fundamental-conjugate representations of sl{M\N') in [30]. 
Open spin chains based on s/(l|2) have been studied in details in e.g. [21,22]. 

The plan of our presentation is as follows: 

We flrst deflne the relevant algebraic objects, R matrix and K matrix, together with their com- 
patibility (Yang-Baxter and reflection) equations and their relevant properties. In particular we 
introduce the various reflection equations which arise in the SP and SNP cases. Let us emphasise 
that all notions introduced in the sl{M) case will be straightforwardly generalised to the sl{M.\N') 
case, albeit with a graded tensor product. 

In a second part we deflne the commuting transfer matrices which can be built in both cases, and 
the local Hamiltonians which can be built from them. Locality requirement leads to considering a 
product of two transfer matrices, resp. soliton-antisoliton and antisoliton-soliton, in the SNP case. 
Once again this construction will be valid, with suitable modiflcations, for the sl^M-lAf) case. 

In a third part we discuss the symmetries of these transfer matrices induced by their respective 
YBE and RE structures, in particular focusing on the connection between the SNP case and twisted 
Yangians. 

We then start the discussion of the analytical Bethe Ansatz formulation for the sl{M) spin chains 
in the SNP case. The derivation of suitable new fusion formulae explicited in Appendix A and B 
makes it possible to get a set of Bethe equations. 

In Section 5 we consider the case of s/(A^|A/') super algebra as underlying algebra. Contrary to 
the previous case it is flrst needed to establish a classiflcation for the reflection matrices based on 
the rational (super Yangian) quantum i?-matrix solution, both for SP and SNP conditions. We then 
establish the Bethe equations for both SP and SNP cases. In the SP case in addition we consider 
spin chains built from general K matrix solutions, in the SNP case we restrict ourselves to diagonal 
K matrices. 

1 Yang— Baxter and reflection equations 

The R and K matrices obey sets of coupled consistency equations together with characteristic prop- 
erties which we now describe. 

1.1 The R matrix 

We will consider in a flrst stage the sl{M) invariant R matrices 

i?12(A) = AI + iPl2 (1.1) 



where V is the permutation operator 

Pl2 = 5^ ^., ® ^,. . (1.2) 

Eij are the elementary matrices with 1 in position (i, j) and elsewhere. 

We define a transposition * which is related to the usual transposition '^ by {A is any matrix): 



A* = V-^ A^ V where { 



' V = antidiag(l, !,...,!), for which 1^2 ^ ^q = 1 

or 

V = antidiagf !,...,!, —1, . . . , — 1 j , for which V"^ = 9o = —1 . 

(1.3) 



M/2 J\f/2 



The second case is forbidden for Af odd. 

This R matrix satisfies the following properties: 

(i) Yang-Baxter equation [1,31-33] 

RuiXl - A2) i?13(Al) i?23(A2) = i?23(A2) i?13(Al) i?12(Al " A2) (1.4) 

(a) Unitarity 

i?i2(A)i?2i(-A)=C(A) (1.5) 

where i?2i(A) = VuRuWVu = Ru'W = ^12 (A). 
(in) Crossing-unitarity 

R{\{X)R%{^\~2tp) =aX + tp) (1.6) 

where p = y and 

C(A) = (A + ^)(-A + t), C(A) = (A + tp){-X + tp). (1.7) 

It obeys 

[A1A2, i?i2(A)] = for any matrix A. (1.8) 

The R matrix can be interpreted physically as a scattering matrix [24, 33, 34] describing the 
interaction between two solitons that carry the fundamental representation of s/(A/'). 

To take into account the existence, in the general case, of anti-solitons carrying the conjugate 
representation of sl{J\f), we shall introduce another scattering matrix, which describes the interaction 
between a soliton and an anti-soliton. This matrix is derived as follows 

R-uiX) = Ri2{X) ■■= R{\{-X-tp) (1.9) 

= R{l{~X-ip) =: Ru{X) = R2i{X) (1.10) 



In the case M = 2 and for 6'o = — 1 {sp{2) case), R is proportional to R, so that there is no genuine 
notion of anti-sohton. This reflects the fact that the fundamental representation of sp{2) = sl{2) is 
self-conjugate. This does not contradict the fact that for A/" = 2 and for 6'o = +1 (so(2) case), there 
exists a notion of soliton and anti-soliton. 

The equality between Ri2{X) and Ri2{X) in (jl.9|) reflects the CP invariance of R, from which one 
also has R12 = Ru, i.e. the scattering matrix of anti-solitons is equal to the scattering matrix of 
solitons. In ()1.9j) . ^i2(A) = (—A — ip)I + iQu, in which Q12 is proportional to a projector onto a 
one-dimensional space. It satisfies 

Q^ = 2pQ and V Q = Q V = OqQ . (1.11) 

The R matrix p.lO|) also obeys 
(i) A Yang-Baxter equation 

RuiXl - A2) Rl3{\l) R23i\2) = i?23(A2) i?13(Al) i?12(Al " A2) (1.12) 



(a) Unitarity 



(Hi) Cms sing -unitarity 



i?i2(A)i?2i(-A)=C(A) (1.13) 



^i2(A) RU-^ - 2^p) = C(A + ip) . (1.14) 



Remark: The crossing-unit arity relation written in the literature usually involves a matrix M = 
V'^V. In our case, M turns out to be 1 for two reasons: (i) the factors q^ of the quantum (trigono- 
metric) case degenerate to 1 in the Yangian (rational) case; (ii) the signs usually involved in the 
super case are in this paper (sectional) taken into account in the definition of the super-transposition 

1.2 The K matrix 

The second basic ingredient to construct the open spin chain is the K matrix. We shall describe in 
what follows two different types of boundary conditions, called soliton preserving (SP) [8-11] and 
soliton non-preserving (SNP) [12-14]. 



1.2.1 Soliton preserving reflection matrices 

In the case of soliton preserving boundary conditions, the matrix ii' is a numerical solution of the 
reflection (boundary Yang-Baxter) equation [3] 

RabiXa - \b) Ka{\a) Rba{K + Afc) Kb{\b) = Mh) Rab{\a + Afe) ir„(A„) RbaiK - \b) , (1.15) 



and it describes the reflection of a soliton on tlie boundary, coming back as a soliton. 

Anotlier reflection equation is required for wliat follows, in particular for the 'fusion' procedure 
described in the appendices 

RabiXa - \b) K-a{\a) i?6a(Aa + Afe) Kb{\b) = Kb{\b) Rah{\a + h) MXa) RbaiK - \b) • (1.16) 

Ka is a solution of the anti-soliton reflection equation obtained from p.l5|) by CP conjugation and 
actually identical to ()1.15|) due to the CP invariance of the i?-matrix. It describes the reflection of 
an anti-soliton on the boundary, coming back as an anti-soliton. 

Equation p.l6|) appears as a criterion for a consistent choice of a couple of solutions Ka and Ka 
of ()1.15j) . yielding the commutation of transfer matrices (hereafter to be defined). 

Graphically, p.l5|) and p.l6|) are represented as follows: 







These K matrices (solutions of the soliton preserving reflection equation p.l5|) ) have been classified 
for sl{Af) Yangians in [35]. This classification can be recovered as a particular case of our proposition 
15.11 where the sl{M. \M) Yangians are studied. Yang-Baxter and reflection equations will indeed take 
the same form albeit with a graded tensor product in the superalgebraic case (see section El for more 
details). Proposition 15.21 then provides the classification of pairs {Ka{\),Ka{\)} which obey (J1.15J) 
and the compatibility equation ()1.16|) . 



1.2.2 Soliton non-preserving reflection matrices 

In the context of soliton non-preserving boundary conditions one has to consider [12-14] the case 
where a soliton reflects back as an anti-soliton. The corresponding reflection equation has the form 



RabiK - K) Ka{\a) Rba{K + \) MXb) = K^,{\) Rab{\a + h) Ka{\a) Rba{K " A;, 



1.17) 



Note that equation ()1.17p is satisfied by the generators of the so-called twisted Yangian [36,37], which 
will be discussed in section 13.21 

Similarly to the previous case, one introduces Kg,-, describing an anti-soliton that reflects back as 
a soliton, satisfying (J1.17p and the consistency condition 

RabiK - \b) K-a{K) Rba{K + A,) Kb{\b) = Kb{\b) Rab{\a + \b) K-a{\a) Rba{K ' \b) ■ (1-18) 

Graphically, (jl.l7|l and p.l8|l are represented as follows: 







The K matrices corresponding to the sohton non-preserving reflection equation are classified in 
proposition 15.31 in the case of sl{}A\M) Yangians, where once again similar Yang-Baxter and refiec- 
tion equations occur. Proposition 15.41 then provides the classification of pairs {i^a(A),i^a(A)} which 
obey ()1.17p and the compatibility equation p.lSp . 

2 The transfer matrix 

We are now in a position to build open spin chains with different boundary conditions from the 
objects K, K, R and R [4]. Our purpose is to determine the spectrum and the symmetries of the 
transfer matrix for the case where soliton non-preserving boundary conditions are implemented. We 
first recall the general settings for the soliton preserving case. 



2.1 Soliton preserving case 

Let us first define the transfer matrix for the well-known boundary conditions, i.e. the soliton 
preserving ones. The starting point is the construction of the monodromy matrices 



r,(A) = T,(A) ir-(A) T„(A), 
T„(A) = T,(A) Kr(A) f,(A) 

and two transfer matrices (soliton-soliton and anti-soliton-anti-soliton) 

t(A) = TiaKiX) %{X), t(A) = TiaK^X) TaiX) 

with 



(2.1) 
(2.2) 

(2.3) 



T,{X)=RaL{X)...Ral{X) 
T-a{X)=RaL{X)...Ral{X) 



f„(A)=i?i,(A)...i?ia(A) 
fa{X)=Ria{X)...RLa{X) 



(2.4) 



The numerical matrices K^ (A), K^ (A) are solutions of p.l5|) . p.l6|) and K^ satisfies a refiection 
equation 'dual' to ()1.15p . 



Rati-Xa + A,) K+{XaY Rbai'Xa - Afc - 2ip) K^ {X,f 

= K+iXbY Rab{-Xa - Ab - 2tp) K+iXaY Rbai-^a + X, 



(2.5) 



The solutions of ()2.5p take the form K^{X) = -K'*(— A — ip), where Ka{\) is a solution of ()1.15|) . 
Similarly ^^^'"(A) satisfies a reflection equation dual to ()1.15p . the solutions of which being of the 
form K^{X) = i^^(— A — ip). In addition, K^ satisfies also a compatibility condition dual to (jl.l6|) . 
Actually, the dual reflection equations happen to be the usual reflection equations after a redefinition 
Ac -^ -Ac - ip- 

From their explicit expression, one can deduce that the monodromy matrices T(A) and T{X) obey 
the following equations: 

Rabija — \b) Ta{\a) Rba{K + Afe) %{\) = %{\) Rabija + Afo) 7^(Aa) Rba{^a — Afc) (2.6) 

Rabija — Afo) Ta{Xa) Rbai^a + Afe) Tb{\b) = T ^{Xb) Rabija + Afe) Ta{\a) Rba{^a — Afc) (2.7) 

RabiK — Af,) Ta{Xa) Rbai^a + Afe) 7^(Afc) = %{Xb) RabiK + Afe) T a{Xa) Rbai^a — Xb) , (2.8) 

which just correspond to the soliton preserving reflection equations p.l5j) and p.lfjj) . One can recog- 
nise, in the relation ()2.6|) . the exchange relation of reflection algebras based on i?-matrix of y{slj\f). 
The form of K^ determines the precise algebraic structure which is involved, see section IHl 

As usual in the framework of spin chain models, the commutativity of the transfer matrices (j2.3p 

[t{X),t{p)] = (2.9) 

[t(A),t(/i)] = [t(A),t»] = (2.10) 

is ensured by the above exchange relations. The first commutator guarantees the integrability of the 
model, whose Hamiltonian is given by 



«-^>' 



(2.11) 

A=0 



the locality being ensured because -R(O) = V. 

The commutators ()2.10|) will be needed so that the fusion procedure be well-defined (see appendices). 

The transfer matrix for K~^{u) = K~{u) = 1 satisfies a crossing-like relation (see for details [14,27]) 

t(\) =t{-X-ip) . (2.12) 

The eigenvalues of the transfer matrices as well as the corresponding Bethe Ansatz equations have 
been derived for diagonal K matrices in [8, 16]. 

2.2 Soliton non-preserving case 

This case was studied in [14] for the sl{3) chain only. Here we generalise the results for any s/(A/'). 
One introduces two monodromy matrices 

UX) = T„(A) ^-(A) t,(A) , 
TaiX) = r,(A) irr(A) t(A) , (2.13) 



and two transfer matrices (anti-soliton-soliton and soliton-anti-soliton) defined by 

t(A) = Tr,i?+(A) Ta{X) , t(A) = Tr,K+(A) r„(A) , (2.14) 

where now 

^a(A)=-Ra2L(A)-Ra2L-l(A) . . .-Ra2(A)-Ral(A) , Ta(A)=-Ri a(A)i?2a(A) • • • -R2L-I a(A)-R2La(A) , 

Ta(A)=-Ra2L(A)-Ra2L-l(A) . . . -Ra2(A)-Ral (A) , Ta(A)=-Ria(A)-R2a(A) . . • -R2L-I a(A)-R2La(A) . 

(2.15) 

Note that, in this case, the number of sites is 2L because we want to build an alternating spin chain, 
which is going to ensure that the Hamiltonian of the model is local. This construction is similar to 
the one introduced in [38] , where however a different notion of R was used. 

The numerical matrices K~ , K^ are solutions of ()1.17j) . p.l8j) . The numerical matrices K'^ and K^ 
are solutions of the following reflection equations: 

Rab{-K + A,) K+{Ky RU-K - At - 2ip) K+{Xt,y 

= K+iXbY Rabi-\a -\b- 2tp) K+iKY Rba{-Xa + Xb) (2.16) 

Rab{-Xa + Xb) K+iXaY Rba{-\a "A, - 2ip) K^ [Xb]' 

= K-+{XhY Rab{-K - A, - 2tp) KtiXaY Rba{~K + \) (2.17) 

Rab{-K + A,) Kt{XaY Rbai-Xa "A, - 2tp) K+ {X^Y 

= K+{X,Y Rabi-Xa - A, - 2tp) K+{XaY Rba^Xa + Afe) (2.18) 

The commutators 

[t(A),t(/i)] = , [i{X),i{p)] = and [t{X),i{p)] = (2.19) 

are ensured by the above exchange relations for K~^ and the relations for the monodromy matrices, 
namely 

RabiK — Aft) Ta{Xa) Rbai^a + Afo) 7^(Ab) = 7^(Ab) Rabija + Afo) Zj,(Aa) RbaiK " Afe) , (2.20) 

Rabija — Afo) Ta{Xa) Rba{Xa + Afo) Tb{Xb) = Tb{Xb) -Rafe(Aa + Afo) T a{Xa) Rba{Xa — Xb) , (2.21) 

Rab{Xa — Xb) Ta{Xa) Rba{Xa + Xb) Tb{Xb) = Tb{Xb) Rab{Xa + Xb) T a{Xa) Rba{Xa — Xb) ■ (2.22) 

The relation ()2.2()|1 has to be compared with the exchange relation of twisted Yangians based on 
i?-matrix of y{slj\f). We will come back to this point in sectional 

In the SNP case, one can show [14] that the transfer matrices for K^{u) = K^{u) = 1 exhibit a 
crossing symmetry, namely 

t(A) = t(-A - ip) , t(A) = t(-A - ip) . (2.23) 

Starting from ()2.15j) and ()2.14j) . the Hamiltonian of the alternating open spin chain is derived as 

[14,17] 

n = -l^t{X) i{X) . (2.24) 

2 aX A=o 



The locality is again ensured by -R(O) = P, and integrability is guaranteed by ()2.19p . 

We can explicitly write the open Hamiltonian in terms of the permutation operator and the -R(O) 
and -R(O) matrices. Let us first introduce some notations: 



d ^ 



dX 



d 



R,, = i?,,(0), R[^ = JTR^A>^) ^_ and R[^ = Vi, tt^.,(A) ^_ • (2.25) 



A=o '' '' dX 



A=0 



After some algebraic manipulations, in particular taking into account that Ttq Vqi Rqi oc I, we obtain 
the following expression for the Hamiltonian ()2.24|) (for a detailed proof see [14]) 

L L-l 

H oc 2_^ i?2j-l 2j -R2J-I 2j + 2_^ R2J+I 2j+2 R2J 2j+2 R2J+I 2j+2 

3=1 i=i 

L-l 
+ 2^ R2J+I 2j+2 R2J-I 2j R2J-I 2j+2 R2J-I 2j+2 -R2J-I 2j R2J+I 2j+2 

i=i 

L-l 
+ 2_^ R2J+I 2j+2 R2J-I 2j R2J-I 2j+2 R2J-I 2j+l -R2J-I 2j+2 ^2]-! 2j R2J+I 2j+2 

i=i 

+ Tro Rq 2L R2L-I 2L Vq 2L-1 Rq 2L-1 R2L-I 2L + Rl2 Ru -^12; (2.26) 

which is indeed local including terms that describe interaction up to four first neighbours. 

It is easily shown, acting on ()2.15p by full transposition, that t(A) oc t^^-^^^{X) provided that 
{K^y oc K^ . Eigenvectors of H in (J2.26J) are determined by sole evaluation of eigenvectors of t{X). 
We shall therefore only need to consider diagonalisation of t{X) in what follows. 

3 Symmetry of the transfer matrix 

In the two (SP and SNP) boundary cases, the use of exchange relations for the monodromy matrices 
allows us to determine the symmetry of the transfer matrix. For simplicity, we fix K^{X) (or K^{X)) 
to be I, leaving K~{X) (or K~{X)) free. 

3.1 Soliton preserving boundary conditions 

In this case, the general form for K~{X) [35] is conjugated (through a constant matrix) to the 
following diagonal matrix: 

K-{X)=i^I + XE with E = diag(+l, ..., +1, -1, ..., -1) . (3.1) 

m n 

In the particular case of diagonal solutions, one recovers the scaling limits of the solutions obtained 
in the "quantum" case in [39]. The monodromy matrix then T(A) generates a B{J\f,n) refiection 
algebra as studied in [40]. Taking the trace in space a of the relation ()2.(ij) . we obtain: 

(A^-A^)[t(A.),T(A,)] = (2tXa-U)[T{Xa),T{X,)] (3.2) 



From the asymptotic behaviour (A -^ oo) of the R matrix 

i?o^(A) = A(I+^Po^) , (3.3) 

we deduce 

(3.4) 



r(A) = A^^+^ f E + ^ (e + E i3o,) + 0(1)") 



where Bqj = Vqj Eq + Eq Vqj. We can write Bqj = ^^ ^^^^ Eap^bj , with i^^^ the elementary matrices 
acting in space 0, and 6° reahzing the gl{m) © gl{n) algebra (in the space j). 
Picking up the coefficient of Ap in the relation ()3.2p . one then concludes that 

:0, (3.5) 

i.e. the transfer matrix commutes with the gl{m) © gl{n) algebra. 



L 

a/3 






3.2 Soliton non-preserving boundary conditions 

We have already mentioned that the monodromy matrix T(A) satisfies one of the defining relations 
()2.2()|) for the twisted Yangians 3^^(A/') [36,37]. The general form for K~{X) is given in proposition 
15.31 it is constant and obeys {K^Y = tK" with e = ±1. 

We need to investigate the asymptotic behaviour (A -^ oo) of the R and R matrices given by 
(jl.ip and (jl.lUp respectively: 

Ro^{\) = a(i + ^Po^) and /?oi(A) = -A(i + ^Poi) , (3.6) 

where we have introduced 

Voi = pl-Qoi- (3.7) 

Accordingly, the monodromy matrices (|2.15j) take the following form 

L 

L 

i=l 
L 



ro(A) = (-AY(i+^E(^o-2. + ^o,2.-i)+0(1) 

fo(A) = (-AY(i+^E(^o.2. + Po.2.-i)+0(1)) (3.^ 



and finally 



where 



_ . 2L 

ro(A) K- %{\) = A^^(ir- + ^ 5^5o. + 0(-)), (3.9) 

So,2i = Vo,2iK- + K~Vo,2i , 5o,2i-l = K-Po,2i-l + Vo,2i-iK" . (3.10) 

10 



Similarly to the previous case, one can show from equation ()2.20j) that (this time for fT =1 only): 



t(Aa),5^5, 



2L 



i=l 



(3.11) 



In this particular case S realises the so{J\f) (or sp{J\f)) generators. When J\f is even, as already 
mentioned, there are two possibilities for the projector Q. More specifically the choice 6'o = 1 in p.3|) 
corresponds to the so{J\f) case, whereas 9o = —1 corresponds to the sp{Af) case. 

Remark: The same construction (open chain with twisted boundary conditions) can be done start- 
ing from the so(A/'), sp{2J\f) and osp{A4\2N') i?-matrix [41]. However, since the fundamental rep- 
resentations of these algebras are self-conjugated, solitons and anti-solitons define the same object. 
Hence, the "twisted" boundary conditions should be equivalent to open chains with "ordinary" 
boundary conditions. Indeed, it has been shown in [41] that boundary reflection equations (defining 
the boundary algebra) and twisted reflection equations (defining the twisted Yangian) are identical. 

4 Spectrum of the transfer matrix 

Our purpose is to determine the spectrum of the transfer matrix for the sl{J\f) case. 

4.1 Treatment of non-diagonal reflection matrices (SP case) 

In the soliton preserving case, the classification of reflection matrices associated to the Yangian 
Y{gl{M)) has been computed in [35]. It can been recovered from the sl{M.\N) case given in propo- 
sition 15.11 Using this classification, it is easy to show the following proposition [42] : 

Proposition 4.1 Let K{X) be any diagonalizahle reflection matrix. D{X), the corresponding diagonal 
reflection matrix, can he written as 

D{X) = U-^K{X)U (4.1) 

where U , the diagonalization matrix, is constant. Let tx(A) = Ti a(Ta{X)Ka{X)Ta{X)) and toi^) = 
TYa(Ta{\)Da{\)Ta{X)) tc thc Corresponding transfer matrices (we set K+ = I). 
Then, ti^(A) and t£){X) have the same eigenvalues, their eigenvectors (say vk and vd respectively) 
being related through 

vk = UiU2...UlVd (4.2) 

Proof: The fact that the diagonalization matrix is a constant (in A) is a consequence of the classifi- 
cation (see [35] and proposition 15. Hi . Using the property ()1.8|1 . one can show that 

tK{\) = U1U2 ...Ul toiX) {U1U2 . . . Ul)-' , (4.3) 

which is enough to end the proof. ■ 

The general treatment (including the super case) for diagonal reflection matrices is done in section 
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I5.4.41 From the above property, this treatment (for diagonal matrices) is enough to obtain the 

spectrum for all the transfer matrices associated to all the reflection matrices (provided they are 

diagonalisable). 

As an illustration of this proposition, we compute the eigenvalues associated to the non-diagonal 

reflection matrix [3, 10] 



/ -X + i^ 

cX + i^ 



KiX) 




V 2kX 











2kX \ 




cX + i^ 
A + ze / 



It is easy to see that K{X) is diagonalized by the constant matrix 





U 



( -\ 

1 



V ^^ 



1 





f \ 





c+1 / 
2c / 



The corresponding diagonal matrix is given by 



with c^ = 4A;^ + 1 



(4.4) 



(4.5) 



L'(A) = cdiag(-A + i^',A + i^',...,A + iO, with ^ 



(4.6) 



in accordance with the classification of refiection matrices. The t/^ (A) eigenvalues as well as the 
Bethe equations are identical to the ones of to (A). They can be deduced from the general treatment 
given in section IB.4.4| taking formally ni = ^2 = A/" = and specifying vti\ = 1. They can also be 
viewed as the scaling limit of quantum groups diagonal solutions [15]. The eigenvectors are related 
using the formula ()4.2|1 . with the explicit form ()4.5j) for U. 



This general procedure can be applied for an arbitrary spin chain, provided the diagonalization 
matrix is independent from the spectral parameter and commutes, see eq. p.8|) . with the i?-matrix 
under consideration. When A'+(A) is not I, this technics can also be used if A'+(A) and K^{X) can 
be diagonalized in the same basis [43]. Note however that the classification does not ensure the full 
generality of such an assumption. 

4.2 Pseudo-vacuum and dressing functions 

We present below the case when soliton non-preserving boundary conditions are implemented with 
the simplest choice K^ = I. Results for more general choices of K~ can be deduced from the 
superalgebra case treated in section IB. 5. 31 
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We first derive the pseudo- vacuum eigenvalue denoted as A°(A), with the pseudo- vacuum being 

2L 

k+) = 6^\+)^ where |+) = 



i=l 







e C^. 



(4.7) 



Vo/ 



Note that the pseudo- vacuum is an exact eigenstate of the transfer matrix ()2.3|) . and A°(A) is given 
by the following expression 

A°(A) = (a(A)6(A))2^(7o(A) + (KA)fe(A))^^ J^ gi{X) + (a(A)6(A))2V-i(A). (4.8) 



1=1 



with 



and 



a(A) = A + z, 6(A) = A, a(A) = a(-A - zp), 6(A) = 6(-A - zp) 



Af-1 



(4.9) 



9iW - 


A + f ' - 


QAf-liX) = 


= 1, for A/" odd 


9iW - 


= gAf-i-ii-X -ip) . 



(4.10) 



M 



We remind that p = y. 

We make at this point the assumption that any eigenvalue of the transfer matrix can be written as 

M-2 

A(A) = (a(A)6(A))2^^o(A)Ao(A) + (6(A)6(A))2^ Y. 9i{X)M>^) + {a{X)h{\)f'^gM-i{mN-i{X) (4.11) 



1=1 



where the so-called "dressing functions" Ai{X) need now to be determined. 
We immediately get from the crossing symmetry ()2.23|) of the transfer matrix: 

Ao(A) = A^_i(-A - ip), Ai{\) = A^_i_,{-X - ip) . 



(4.12) 



Moreover, we obtain from the fusion relation ()A.10|) the following identity, by a comparison of the 
forms 1)4.111) for the initial and fused auxiliary spaces: 



Ao{\ + ip)A^^i{\) = l. 
Gathering the above two equations ()4.12j) . ()4.13|) we conclude 

Ao{\)Aoi-X) = l. 



(4.13) 



(4.14) 



Finally from equations ()B.13|) important relations between the dressing functions are entailed for 
both soliton preserving and soliton non-preserving boundary conditions. In particular. 



Af-l 



Yl Ai{X + i{Af - 1) - il) = 1. 



(4.15) 
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Taking into account the constraints ()4.12j) . ()4.14|) and ()4.15|) one derives the dressing functions: 

MX) = U ' ~' ' "' 



(i)^i \_ \(i) 



--1 X + Xj + J X- Xj +1 

y A + Af + f A-Af + f 

X TT ^y— ? ^ ly— ? ^, l</< 4.16 

M A + Arm + t A-Afm + t' - 2 ^ ^ 

together with the property Ai{X) = y4_^_i_;(— A — ip), and, for A/" = 2r2 + 1: 

AJX) = n \ , -^^ 7-1 fi -, 4.17 

f}, X + xf^ + f X-Xf + f A + Af) + f + 1 A-Af) + f + 1 

Note that the dressing does not depend on the value of ^o- 

The numbers M^''^ in ()4.17|) are related as customary to the eigenvalues of diagonal generators Si 
of the underlying symmetry algebra (determined in the previous section), namely 

Si = ^M(o)-M«, ^; = M('-i) - MW with Si = ^iEu-Eji), 1 < / < :^^ (4.18) 

with M(°) = 2L. 

Recall that for the sl{J\f) case the corresponding numbers M^^\ see e.g. [11], are given by the following 

expressions 

Ell = M^^-^^ - M^^\ / = 1,...,AA (4.19) 

with M^^^ = 2L and M*^''*^-' = 0. If we now impose M*^') = M'^^~^'^ and consider the differences 
Ell — Eji, we end up with relations (|4.18|) . in accordance with the folding of sl{J\f) leading to so{Af) 
and sp{Af) algebras. 

4.3 Bethe Ansatz equations 

From analyticity requirements one obtains the Bethe Ansatz equations which read as: 
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4.3.1 sl(2n + 1) algebra 

mW Af(2) 



1 = - n -^(^? - 4') ^^(^? + ^f ) n n -i(^? - ^r^) -^(^i' + 4 



j = l r=±l j=l 



l = 2,...,n-l, 

i=i 

X n e_i(Af)-Af-'V-i(A?^ + Af"'^ (4.20) 

i=i 

where we have introduced 

A -I- — 

It is interesting to note that equations ()4.2()j) are exactly the Bethe Ansatz equations of the osp(l|A/'— 
1) case (see e.g. [14,19]). 

4.3.2 sl(2n) algebra 

e.{^^Y = - n ^^(^' - ^') ^-^ + ^r) n ^-^^ - ^f ) ^-^^ + ^f ) ' 

3=1 i=i 

l = -ne,(Af-Af)e,(A,'" + A;") J] 11 e-,(A<" - A?-') e_.(A.'" + A-'-') 

i = l T = ±l j = l 

/ = 2,...,n-l, 



-,..., .1/ -1-, 

j 



e-.o(ArO = - n ^2(^S"^ - ^i ) e2(Ar^ + Aj. 



X 



The Bethe Ansatz equations are essentially the same as the ones obtained from the folding of the 
usual sl{M) Bethe equations (see e.g. [8,11]) for M^'^ = M'^-'^~'\ It can be realised from the study of 
the underlying symmetry of the model that this folding has algebraic origins, as mentioned previously. 
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5 sl{A4\Af) superalgebra 

In this section, we generalise the previous approach on (SP and SNP) boundary conditions to the 
Z2-graded case based on the sl{A4\Af) superalgebra. 

5.1 Notations 

The Z2-gradation, depending on a sign 6*0 = ±, is defined to be (—1)'*^ = Oq for i an sl{Ai) index 
and (—1)'*' = —00 an sl{J\f) index. 
The sl{Ai\f/) invariant R matrix reads 

RuiX) = XI + iPi2 , (5.1) 

where P is from now on the super-permutation operator (i.e. X21 = PX12P) such that 

M+M 

P= 5^(-l)[-''l^.,®% (5.2) 

The usual super-transposition -^ is defined for any matrix A = ^. A"^^ Eij, by 

A'^ = J](-i)Wb1+b1 A^^ E,, = Y, {^1^' E,, . (5.3) 



t 



As for the sl{J\f) case, we will use a super-transposition * of the form 

A' = V~^A^V. (5.4) 

The convention for 9q and the expression of V are chosen accordingly to the selected Dynkin diagram. 
Let us recall that for a basic Lie superalgebra, unlike the Lie algebraic case, there exist in general 
many inequivalent simple root systems (i.e. that are not related by a usual Weyl transformation), 
and hence many inequivalent Dynkin diagrams. This situation occurs when a simple root system 
contains at least one isotropic fermionic root. For each basic Lie superalgebra, there is a particular 
Dynkin diagram which can be considered as canonical: it contains exactly one fermionic root. Such 
a Dynkin diagram is called distinguished. In the case of s/(A1|A/'), it has the following form: 

(y^ ® (>;0 

M-1 M-1 

In the case of sl{Ai\2n) superalgebras, there exists a symmetric Dynkin diagram with two isotropic 
fermionic simple roots in positions n and Ai + n: 

c>^jO — ® — Q::sS^ — ® — ^^:u^ 

n— 1 M — l n—1 
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The generalization of the Weyl group for a basic Lie superalgebra gives a method for constructing 
all the inequivalent simple root systems and hence all the inequivalent Dynkin diagrams. For more 
details, see [44-47]. 

(i) Distinguished Dynkin diagram basis 

In this case, we consider that the sl{Ai) part occupies the 'upper' part of the matrices and 
corresponds to bosonic degrees of freedom, whereas the sl{Af) part occupies the 'lower' part 
and corresponds to fermionic degrees. More precisely, the gradation takes the form: 

-1)^^^ = ^ ' for 1<.<A^ 

^ ^ -1 for M + l<i<M+X, 



and the matrix V reads 



y_ , Vm 



° ' (5.6) 



vv 





In the above formula, Vm (resp. VV) is the Ai x A4 (resp. A/" x A/") matrix given in ()1.H|1 for 
00 = +1. 

(ii) Symmetric Dynkin diagram basis 

As in the sl{Af) case, one has to take Al or A/" even. Note that for the odd-odd case no 
symmetric Dynkin diagram exists and consequently no twisted super- Yangian [48]. Here, we 
choose jV to be even: A/" = 2n. The sl{Ai) part lies in the 'middle' part of the matrices and 
corresponds to fermionic degrees of freedom, whereas the sl{Af) part occupies the 'upper' and 
'lower' part and is associated to bosonic degrees of freedom. Correlatively, 6*0 = —1 in this 
case. The gradation is given by 

[i] _ J 1 for l<i<n and M + n + 1 < i < M + Af 

[-1 for n+l<i<M + n, 

while 



V = antidiagf 1, . . . , 1 , -1, . . . , -1 j . (5.8) 

n+M n 

We will mostly use the distinguished Dynkin diagram basis in the soliton preserving case, and the 
symmetric one in the soliton non-preserving case. In both cases, the i?-matrix obeys the properties 
stated in section IlT| with -Ri2(A) = R%{—X — ip) and 2p = Oq{M. — A/"). The i^-matrices will obey 
the defining relations stated in section 11.21 and the properties of the transfer matrix (section |2)) also 
hold; the tensor product is now Z2-graded. 

5.2 Classification of reflection matrices for y{J^\Af) 

This section is devoted to the classification of reflection matrices for the super- Yangian y{A4\N') 
based on s/(A/l|A/'), both for soliton preserving (prop. 15.11 and l5.2|) and for soliton non-preserving 
boundary conditions (prop. l5.Hl and |5.4j) . 
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5.2.1 Soliton preserving reflection 

Proposition 5.1 Any bosonic invertible solution of the soliton preserving reflection equation (RE) 
RuiXl - X2)Ki{X,)Ru{Xl + \2)K2{\2) = K2{\2)Rl2{\l + A2)iri(Ai)i?i2(Ai - A2) 

where Ri2{\) = \l + i P12 is the super- Yangian R-matrix, takes the form K{X) = U (?^ I + A E) U~^ 
where either 

(i) E is diagonal and E^ = I (diagonalisahle solutions) 
(a) E is strictly triangular and E^ = (non- diagonalisahle solutions) 

and U is an element of the group GL{Ai) x GL{N'). The classification is done up to multiplication 
by a function of the spectral parameter. 

Proof: Firstly it is obvious tliat for any solution K{\) to the RE, and for any function /(A), 
the product f{\)K{\) is also a solution to the RE, so that the classification will be done up to 
multiplication by a function of A. 

Expanding the reflection equation, one rewrites it as: 

[K^iX,), K2{X2)] = 2(Ai + A2)(i^2(Al)i^l(A2)-i^2(A2)i^l(Al)) 

+2(Ai - A2) (iri(Ai) K,{X2) - K2{X2) K2{Xi)) (5.9) 

One then considers the RE with Ai and A2 exchanged, and sums these two. After multiplication by 
P12, one gets (for Ai 7^ A2): 

[i^i(Ai), i^i(A2)] = -[i^2(Ai), i^2(A2)] (5.10) 

the only solution of which is [K{Xi), K{X2)] = 0. In other words, the matrices K{X) at different 
values of A's are diagonalisable (or triangularisable) in the same basis and must satisfy 

(Ai + A2) {K2{X,) Ki{X2) - K2{X2) K,{Xi)) + (Ai - A2) (i^i(Ai) i^i(A2) - K2{X2) K2{X,)) = (5.11) 

By setting A2 = — Ai, we get K{X)K{—X) = k{X) I for some function A;(A). If one now considers the 
case of invertible matrices, and since we are looking for solutions up to a multiplicative function, we 
can take K{X)K{—X) = I, a condition which is generally assumed for reflection matrices. 

We first consider the case where these matrices can be diagonalised: K{X) = U D{X) U^^, where 
f/ is a group element of GL{m) x GL{n). Projecting the RE on the basis element En ^ Ejj, one gets 

(Ai + A2) (rf,(Ai) rf,(A2) - d,{X2) rf.(Ai)) = (Ai - A2) (rf,(A2) rf,(Ai) - rf,(Ai) rf,(A2)) (5.12) 

where -D(A) = diag((ii(A), d2{X), . . . , dm+n{X)). Since K{X) is supposed invertible, all the d/s are not 
zero, and we consider 

,«(A) ^ 1^ (5.3) 
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which obeys 

(x + y) \^q{y) - q{x)j + {x - y) [q{x)q{y) - 1 j = . (5.14) 

The solution to this equation is q{x) = — fr^ where ^ is some complex parameter (including ^ = oo), 
so that, considering gji(A), we get 



dj{X) = --—fd,{X), Vj (5.15) 

Requiring qij{X) to obey the equation ()5.14|) shows that one must have 

dj{\) = e^\ + i^ with ej = ±l, Vj (5.16) 

where we have used the invariance under multiplication by a function. This yields the form (i), with 
E = diag(ei,...,em+„). 

We now turn to the case K{X) = UT{X) U^^ where T(A) is triangular. The projection of the RE 
on Eii ® Ejj shows that the diagonal part of T(A) is still of the form {i). We distinguish two cases: 
E has two different eigenvalues (which are ±1), or E is proportional to I (and then the diagonal of 
T(A) is also proportional to I). 

If E has two eigenvalues, we project, in the first auxiliary space, on two diagonal elements Ejj 
and Ekk associated to these eigenvalues: 

(Ai + As) {{ti ± A2)T(Ai) - (z^ ± Ai)T(A2)) = (Ai - A2) (t(Ai)T(A2) - {ti ± \iM ± A2)) (5.17) 

The difference and the sum of these equations read: 

A2T(Ai)-Air(A2) = A2-Ai, Ai + A2 ^ (5.18) 

^^(Ai + A2) (t(Ai) - T(A2)) = (Ai - A2) (T(Ai)r(A2) + C' - A1A2) (5.19) 

From equation ()5.18p . one gets 

r = r = To, I.e. T(A)=I+ATo (5.20) 

Ai A2 

where Tq is a triangular matrix. Plugging this solution in equation ()5.19|) . we obtain 

[i^ - l)(Ai + A2) To = AiA2(T2 - I) + [e + 1) I, VAi, A2 (5.21) 

whose only (constant) solution is of the form {%) with i^ = 1. 

We are thus left with the case where the diagonal of T(A) is proportional to the identity matrix: 
T(A) = I + >S'(A) with 5'(A) strictly triangular. Projecting once more on a diagonal element in the 
first auxiliary space, we obtain 

2 [X2 5(Ai) - Ai 5(A2)) = (Ai - A2)5(Ai)5(A2) (5.22) 

Ai(2I + 5(A0) A2(2I + 5(A2)) ^ ' ^ 
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where a is strictly triangular. We therefore have T(A) = I + 2A(t(I — \a)~^. With this form for 
T(A), the RE rewrites (cxi — cr2)cria2 = 0, whose solution (for strictly triangular matrices) is given by 
cr^ = 0. Using this property, we get the solution {ii). ■ 

Note that the solutions given in this proposition are all of the form K[X) = i^I + XS with S"^ = 1 
or ^2 ^ (taking £ = UEU'^). 

Proposition 5.2 Given a solution K{X) = iC,I + XS to the soliton preserving RE 

RuiXi - A2) i^i(Ai) i?2i(Ai + A2) K2{X2) = K2(A2) RuiXi + A2) iri(Ai) i?2i(Ai - A2) , (5.24) 

and a solution K{X) = i^' I + X S' to the anti-soliton preserving RE identical to \5.24^ , the compati- 
bility condition 

RuiXi - A2) KiiXi) i?2i(Ai + A2) i^2(A2) = i^2(A2) RuiXi + A2) i^i(Ai) i?2i(Ai - A2) (5.25) 

zs solved by £' = S' and ^ + C = Oq^^ StiS. 

Proof: Straightforwardly, equation ()5.25|) is equivalent to 

^1^12^2 = S2Qi2S'l (5.26) 

2{^ + 0[^2,Ql2] = [S2,Ql2S2Ql2] (5.27) 

The first equation yields £' = £^. Using (5i2^2<5i2 = ^0 *^2 Q12 Str£^ one gets the relation between 
C and f . ■ 

5.2.2 Soliton non-preserving reflection 

Proposition 5.3 Any bosonic invertible solution of the soliton non-preserving RE 

RuiXi - A2) K^{X^) i?2i(Ai + A2) i?2(A2) = i?2(A2) RuiXi + A2) K^{X^) i?2i(Ai - A2) (5.28) 

where Ri2{X) = AI + i P12 is the super- Yangian R-matrix, is a constant (up to a multiplication by a 
scalar function) matrix such that K^ = ±K. 

Proof: Writing the R and R matrices in terms of I, Pu and Q12, and taking the part of ()5.28j) which 
is symmetric in the exchange of Ai and A2, yields the following equation 

Ki{Xi)QuKi{X2) + K,{X2)Qi2Ki{Xi) = K2{Xi)Qi2K2{X2) + i^2(A2)gi2i?2(Ai) (5.29) 

In the same way, exchanging the role of spaces 1 and 2 from the original equation, one gets 

i^i(Ai)Qi2i^2(A2) + i^2(Ai)Qi2i^i(A2) = i?2(A2)Qi2i^i(Ai) + K^{X2)QuK2{X,) (5.30) 

Transposing both equations (j5.29|) and (|5.3(J|) in space 1 and eliminating P12, one gets after some 
algebra 

K\X2) = /(Ai, A2)/^(Ai) V Ai, A2 (5.31) 

from which the final result follows. ■ 
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Proposition 5.4 Given a solution Ki to the soliton non-preserving RE 

i?i2(Ai - A2) i?i(Ai) i?2i(Ai + A2) K2{X2) = K2{X2) i?i2(Ai + A2) i^i(Ai) i?2i(Ai - A2) , (5.32) 
and a solution Ki to the CP-conjugate RE identical to l\5. ,?g)) . the compatibility condition 

RuiXi - A2) i?i(Ai) i?2i(Ai + A2) K2{X2) = K2{\2) RuiXi + A2) KiiXi) /?2i(Ai - A2) (5.33) 
is solved by Ki oc (Kij 
Proof: Straightforward. ■ 

5.3 Pseudo- vacuum and its dressing 

We can determine explicitly the eigenvalue Ao(A) of the transfer matrix (defined as in section |21) 
acting on the pseudo- vacuum | ti;+) (which is always bosonic in our conventions). We take here 
K^ = I (resp. K^ = I), whilst cases with non trivial K^ (resp. K^) are studied in section 15.4.41 
(resp. I5.5.3|l . Ao(A) is given by the following expression 

M+Af-2 

A°(A) = «(A)^(7o(A) + /3(A)^ Yl (-1)''^"^KA) + 7(A)^(-l)[-^+-^^A4+Ar-i(A) (5.34) 

1=1 

where for (using the notation given in ()4.9|) ): 

(i) Soliton preserving boundary conditions with L sites (distinguished Dynkin diagram) 

a{X) = a\X), (3{X) = 7(A) = ^^(A) (5.35) 



and 



A(A + ^^^^) 

gi(X) = ^ .„_/,, , l = 0,...,M-l 

^ ' (A + f)(A + ^)' 



2'^ I 2 

i{M-M)^ 



^KA)= ^^ , .(2^-^-1) ..; , .Wo^ ' l = M,...M+U-l (5.36) 

{A H 2 )\^ "^ 2 ) 



(ii) Soliton non— preserving boundary conditions with 2L sites (symmetric Dynkin diagram) 

2 /_\2 / \2 



a(A) = (a(A)6(A)) , /3(A) = (6(A)6(A)) , 7(A) = (a(A)6(A)) (5.37) 



and 



A + |(p + l) M+M-l 

9AX) = —^z ' < / < 



A + f - 2 

g M+^/-l (X) = 1, ii M+Af odd 

9iiX) = gx+M-i-ii-X-ip). (5.38) 

We remind that Oq = —1 in that case. 
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From the exact expression for the pseudo-vacuum eigenvalue, we introduce the following assumption 
for the structure of the general eigenvalues: 

A(A) = a{\fgo{\)A,{\)+P{\f Y. {-^f^'^9iWM\) 

1=1 
+ 7(A)^(-l)[-^+^-il^^+^_i(A)A^+^_i(A) (5.39) 

where the dressing functions Ai{X) need to be determined. The basic constraints that they have to 
satisfy are the fusion and crossing equations as well as analyticity requirements. 

5.4 sl{A4\J\f) with soliton preserving boundary conditions 

From the analyticity of A(A), one gets 

^K-^) = A-ii-'-^), l = l,...,M-l, 
A2M-i{-^) = A2M-i-i{-^). l = M-M+l,...,M-l (5.40) 

Gathering the constraints (j2.12p . ()A.10|) and (J5.4UJ) . one can determine the dressing functions, i.e. 

MX) = n ' "' ' "' 



J=l -^ + \' + 2 -^ ~ \- +2 



AlW = n 



*^" A + Af + f+z A-Af + f + z 
j—i A + A • ~\~ 2 A — A • ~\~ 2 



i=l ^ + \- + 2 + 2 ^ ~ \- ^ 2 



) 

^KA) = n 



*^*' A + Af +2A^-f-^ A-Af +z.M-f-^ 
,=1 A + Af + 2A^ - f A - Xf +iM-'^ 

/ = A<,...,A^+A/'-l (5.41) 
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5.4.1 Bethe Ansatz equations for the distinguished Dynkin diagram 

From analyticity requirements one obtains the Bethe Ansatz equations, 

Af(i) Af(2) 

3=1 J=l 

1 = - n ^^^ - 4') -^(^? + >?) n n -i(^? - ^r^) -^(^j' + m'"^^) 

3=1 r=±l i=l 

/ = 2,...,Al-l,A< + l,...,Al+A/'-2 

X n eilAS^^-Af^^VilAl'^^ + Ar^^^) 
1 = - n e_,(AS^+^-^) - ;^(^4-Ar-2)^ e_,(Al^+^-^) + Af +^-^)) 

X n e2(A(^+-^-^) - Af +^-^)) e2(A;^+-^-^) + Af +-^-^)) (5.42) 

We recover here for A^ = 2 and A/" = 1 the Bethe Ansatz equation of the supersymmetric t — J 
model which corresponds to the s/(2|l) case [49]. 

5.4.2 Bethe Ansatz equations for arbitrary Dynkin diagrams 

We wrote above only the dressing functions that correspond to the distinguished Dynkin diagram. It 
is however possible to construct the g and dressing functions for all the inequivalent Dynkin diagrams 
of s/(Al|A'). 

The inequivalent Dynkin diagrams of the s/(Al|A^) superalgebras contain only bosonic root of 
same square length ("white dots"), usually normalized to 2, and isotropic fermionic roots ("grey 
dots"). A given diagram is completely characterized by the p-uple of integers Q < rii < . . . < rip < 
M. + M labelling the positions of the grey dots of the diagram. Formally, we define rio = and 
np+i = Al + A/" although there is actually no root at these positions. Such a diagram defined by the 
p-uple (nj)j=i p corresponds to the superalgebra s/(Al|A/') with 

M. = y, ^i ~ /, ^i ^^d M = 2, ^i ~ /_, ^i ■ (5.43) 

i odd i even i even i odd 

i<p+l i<p+l i<p-\-l i<p-\-l 

The g functions have a form similar to (j5.36|) . with a change of increasing or decreasing behaviour 
of the poles each time a grey (fermionic) root is met. Indeed 

A fA + HMzIfi] 

^'^^) = 7Y— WY-TTT^T^ ' ^ = 0,...,A1+A--1 (5.44) 

(A + |(5/) (A + f(5i + l)) 
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where ^o = whilst the 6i ioi I = 1, . . . , A4 + JV — 1 are obtained by iteration 

6i-i if I = Ui for some i 

^i = { 5i-i + 1 if n2i < / < n2i+i for some i (5.45) 

5i^i — 1 if n2i^i < I < n2i for some i 

The Bethe Ansatz equations read, for £=1,...,A^+A/'— 1 and i = 1, . . . , M *^^) 

M+^f-l M(*) |'g^Q(l)\2L c=i 

(i-K«.)) n n^("^-"^>(^s'^-4'^)^("^-"^>(^s'^+^r)= ^ (5.46) 

where (a^, at) is the scalar product of the simple roots associated to the chosen Dynkin diagram. 

5.4.3 Bethe Ansatz equations for the symmetric Dynkin diagram 

We give the useful example of the symmetric Dynkin diagram for which A/" is even, with the indices 
ordered as in ()5.7p . The g functions are in this case 



9iW = —, v,/i, , / = 0,...,Ar/2-l 

(A + f)(A + ^) 

9iW = .,,, , ,, ' — TTTT^, l = Ar/2,...,M+Ar/2-l 



9iW = 7r-W^mV^±wkM±W-, ^ l = M+Ar/2,...,M+Ar-l (5.47) 

and it is straightforward to get the Aj's. The Bethe Ansatz equations take the form: 

ei(A«)- = - n e2(A?^ - A«) e,(Af) + A«) J] e-i(A!^) - Af ) e_,(A« + Xf) , 

1 = - n ^^(^? - 4') ^^(^? + ^f ) n n -i(^? - ^r^) -^(^j' + ^"^') 

j=l T=±l i=l 

l = 2,...,M+Ar-2, /^^, :^ + A< 

1= n e.(Af-Ar))e,(Af + Ar)) J] e_,(A« - A^^^) e_,(Af) + A^^)) 

, AT AT ,^ 
2 ' 2 

1 = - n e_i(A;^+-^-^) - Af +^-^)) e_i(A(-^+^-^) + Af +-^-^)) 

X n e2(AS^+-^-^) - Af +^-^)) e2(Af +-^-^) + Af +-^-^)) (5.48) 

i=i 
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The indices of the e's in the products are the entries of the Cartan matrix corresponding to the 
chosen Dynkin diagram, in accordance with the results obtained in the closed chain case [20,29,50]. 

5.4.4 Non trivial soliton preserving boundary conditions 

We come back to the distinguished Dynkin diagram basis and implement non trivial soliton preserving 
boundary conditions. From the classification given in section 15.2.21 we know that K~{X) is always 
conjugated (by a constant matrix U) to a diagonal matrix of the form 



K{\) = diag(a, . . . ,a, (3, . . . , (3, (3, . . . , P,a, . . . ,a) (5.49) 




As in the section 14.11 it is easy to see that the spectrum and the symmetry of the model depend 
only on the diagonal ()5.49p . and not on U. Indeed, when considering two reflection matrices related 
by a constant conjugation, the corresponding transfer matrices are also conjugated. Thus, property 
15.11 ensures that it is enough to consider diagonal K{X) matrices to get the general case. Such a 
property, which relies on the form of the i?-matrix, is a priori valid only in the rational sl{J\f) and 
sl{M\J\f) cases. 

For a diagonal solution ()5.49|1 with mi + 1712 = Ai, ni+n2 = Af, a{X) = —A + i^, j3{X) = X + iC,, and 
the free boundary parameter ^, one can compute the new form gi{X) of the ^f-functions entering the 
expression of Ao(A), the new pseudo- vacuum eigenvalue. They take the form: 

9iW = i-X + t09iiX), / = 0,...,mi-l 

gi{X) = {X + i^ + imi) gi{X), I = mi, . . . ,M + n2 - I 

gi{X) = {-X + i^-im2 + in2)giiX), l = M+n2,...,M+^f-l (5.50) 

where gi{X) are given by ()5.H(Jj) . The dressing functions ()5.41|1 keep the same form, but the Bethe 
Ansatz equations ()5.42|) are modified (by K~{X)), so that the value of the eigenvalues A(A) are 
different from the ones obtained when K{X) = I. 

The modifications induced on Bethe Ansatz equations ()5.42j) are the following: 

- The factor — e^/^„^(A) appears in the LHS of the mi*'' Bethe equation. 

- The factor —e2£+mi-m2-n2W appears in the LHS of the {M. + 112)*'^ Bethe equation. 

5.5 sl{M.\J\f) with soliton non-preserving boundary conditions 

From equations of the type ()B.13|1 for the supersymmetric case relations between the dressing func- 
tions are entailed for both soliton preserving and soliton non-preserving boundary conditions. In 
particular, for the case that corresponds to the symmetric Dynkin diagram one obtains (A/" = 2n, 
while A4 can be even or odd), 

n— 1 n— 1 M — 1 

J] MX - il) J] A^+„+/(A + iM- i{n - 1) + ^/) = J] Ar,+i{X - i{n - 1) + il). (5.51) 

1=0 l=G 1=0 

In fact the latter equation is only necessary for the soliton non-preserving boundary conditions. 
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5.5.1 Dressing functions 

As already mentioned we consider here the R matrix, that corresponds to the symmetric Dynkin 
diagram. Note that the sl{2n\/iA) case is isomorphic to sl{JiA\2n) and entails the same dressing 
functions and BAE. 

From the constraints (J4.12J) . (j4.14|) . (jS.Slj) . we conclude that the dressing functions take the form: 

MW \ , \(1) i \ \(1) 



^o(A) = n 



A + A . ' - ^ A - A . ' - 



2^ 

•j 2 " "j 2 



j=i A + Aj. + f A - A^- + I 
^ ' |i A + Af + f A-Af + f 

i=l ^ + ^i + 2 + 2 ^ ~ ^i +2+2 

M(') \ , \ (0 , ■ il ■ \ \(l) . ■ il 

^KA) = n 



,=1 A + \f + zn - f A - Af + zn - f 
^^'"^'A + Ar^) + zn-| + ^ A-A^^ + zn-f + l , ^ 



X 



n ^ ^ -^ 2 ' 2 :: 2^ _iz 2 ' 2 ^</<^, 
,., A + Af^m-f-t A-Af^m-f-t' ^-^<^ + 



and Ai{X) = A_M+2n-i-i{—^ — ip), and for A4 = 2m + 1 

Af (*' \ , \ (fe) , ■ ik ■ \ \(k) , ■ ik 

TT- A + A . + zn - f - 2 A - A • + zn - f 
^fe(A) =11 



,=1 A + Af ^ + in - f A - Af ^ + in - f 



I 



A + Af ^ + ^n - f + i A - \f +in-i^+ , ^ ^ 

X -OA 7^^ -, k = m + n 5.53 

A + Af + zn-f-i A-Af +m-f-|' 

Recall that the sl{2m + l|2n + 1) case is not examined because there is no symmetric Dynkin 
diagram and consequently the 'folding' of the algebra can not be implemented. In any case it is 
known [48] that the twisted super- Yangian does not exist for sl{2m + l|2n + 1). 

5.5.2 Bethe Ansatz equations 

From the analyticity requirements one obtains the Bethe Ansatz equations which read as: 
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A. sl(2m + l|2n) superalgebra 

M(i) Af(2) 

e,(A:")- = - n -^(Al" - A?') -^(A'" + Af ) n -.(^l" - Af ) e_,(Af' + Af ) . 
i=i i=i 

j = l T = ±l J = l 

/ = 2, ...,?T, + m — 1, l^n 
1= J] ei(Af)-Af+'))ei(A;") + Af+')) J] e_i(Af ^ - Af"')) e_i(A(") + Af"')) 

i=i 

X J] e„i(Af)-Af"^))e_i(Af^ + Af-^)), k = m + n (5.54) 

i=i 

Note that the equations (|5.54p are the Bethe Ansatz equations of the osp{2m + l|2n) case (see 
e.g. [19]) apart from the last equation. 

B. sl(2m|2n) superalgebra 

The first n + m — 1 equations are the same as in the previous case see equation (j4.2(jp . but the last 
equation is modified, with again k = m + n, to 

ei(Af )) = - n ^^(Af ^ - Af) ^2(Af ) + Af ) n ^-i(\^'^ - ^r') ^^i^^ + ^t'') ■ (5-55) 

Notice that the Bethe Ansatz equations as in the non supersymmetric case, are essentially the ones 
obtained from the folding of the symmetric sl{A4/J\f) Bethe equations (j5.48|) for M*^'^ = M^-^~^-'''^^~^\ 
This folding has algebraic origins as can be realised from the study of the underlying symmetry of the 
model (see ()4.18|) . ()4.19p ). Indeed only half of the sl{A4\J\f) generators survive after we impose the 
soliton non-preserving boundary conditions, and these are exactly the generators of the osp{Ai\Af) 
algebra. 

5.5.3 Non trivial soliton non-preserving boundary conditions 

We generalize the above approach to the diagonal case with e = 1 of the classification given in 
proposition 15.31 

K'{\) = diag(/i;i, . . . , kM+j\f) with kM+M+i-j = kj . (5.56) 

We will consider only invertible K'^ matrices, so that ki ^ \/i. 
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The (^-functions entering the new pseudo-vacuum eigenvalue are modified in the following way: 

giiX) = h+,giiX), 0<l<^i±^^ (5.57) 

where gi{X) are given by ()5.H8j) . The remaining g are defined by requiring the crossing relation 

9M+Af-ii->^ - W) = 9iW ■ (5.58) 

The dressing functions (J5.52J) and (|5.53|) keep the same form, but the LHS of £*'^ Bethe Ansatz 
equation (given in ()5.54p and ()5.55|) ) is multiplied by fc^/fc^+i. 

Acknowledgements: This work is supported by the TMR Network 'EUCLID. Integrable models 
and applications: from strings to condensed matter', contract number HPRN-CT-2002-00325. 

A Survey of fusion 

We present here a brief review on fusion for systems with boundaries [51]. In particular we present 
the fusion procedure for open systems without crossing symmetry developed in [14,16]. 
To cover both the SP and SNP cases, we define the action * by: 

R:, = Rab, Kb = Rah. K: = Ka, K*, = K^, for SP b.c. 

R:, = Rab, Kb = Rah, K = Ka, K = Ka, for SNP b.C. (A.l) 

The K* and K^ matrices are solutions of the refiection (boundary Yang-Baxter) equation [3] 

RMa - h) KiK) RISK + h) Ki{\) = Ki{\) Ri,{K + \) K{K) RUK - K), (A.2) 

and they are related by the constraint 

RabiXa - Xb) Kl{Xa) Rl{Xa + Afe) K^iX^) = K^iX,) R*Ma + Afc) K^iXa) KaiXa " Afe) . (A.3) 

These equations unify the equations p.l5|) - ()1.16|) with p.l7|) - ()1.18|) . 

The starting point for both cases (SP and SNP) is the observation that the ^-matrix ()1.9jl 
at the special value A = —ip yields a one-dimensional projector onto the one-dimensional s/(A/')- 
representation present in the decomposition TV ® A/" = 1 © (A/"^ ~ !)• 

P^h = ^ Kh = ^Qah = p;, . (A.4) 

This is related to the fact that the R matrix describes the scattering between soliton and anti- 
soliton. Accordingly, the (A/"^ — l)-dimensional projector is 

Pth=^- P-ah- (A.5) 

Note that the soliton-soliton R matrix at A = — i provides a projector onto a A/"- dimensional space, 
reflected in A/" ® A/" = A/" © (A/"^ —A/"). Thus, one needs the R matrix, even in the SP case, hence the 
introduction of t(A) in both cases. 
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We will formulate the fusion procedure for both types of boundary conditions we mentioned. 
We introduce the fused i?-matrices 

R<-ab>iW = Pn, KiW ^w(A + ^P) PL R<b-a>iW = Kn KiW Kii>^ + ^P) PL (A.6) 

and 

^t<..>(A) = P^, RU\ - ip) RlaW P-L R*i<b-a>W = Pm ^L(A - tp) R^X) PL (A.7) 

They satisfy generalised Yang-Baxter equations with fused indices. Similarly, we use the reflection 
equation ()A.2|1 and its dual to obtain the fused K matrices 



K~,i.>W = Pt,Kl-{\)RU2\ + ip)Kl{X + ip)P, 



+ 



^S.>(A) = nt^r(A)i?L(-2A-3^p)ir*+(A + zp)P+. (A.8) 

Both fused matrices obey generalised reflection equations of the type p.l6|) and its 'dual' (for more 
details we refer the reader to [16,51]). In an analogous way we obtain the K*t matrices by fusing 
the spaces a and b. Now that the fused R and K* matrices are available, we operate the fusion of 
the transfer matrix ()2.3p . The fused transfer matrix is defined by 

MX) = Tr,, {i^*+,>(A) T<,,>(A) ir*-,>(A) f *,,>(A + ^p)} , (A.9) 

where the fused T matrices are obtained using ()2.4|1 with fused R matrices ()A.6j) and ()A.7j) . After 
some algebra (see e.g. [51]) we end up with: 

tpiX) = C(2A + 2zp) t(A) t(A + tp) - A[K*+{X)] 5[T{X)] A[K*-{X)] 6[f*{X)] , (A.IO) 

where C* = C ^iid C* = C i^^ the SP case, while C* = C ^ind C* = C iii the SNP case. Furthermore the 
'quantum determinants' are (when we fuse the spaces a and b) 

5[TiX)] = Tr,5{P^,T,(A)T,(A + ^p)} 

6[f*iX)] = TT,,{P-~,f:iX)f*iX + zp)] 

A[K*-iX)] = Tr,, |P,-, ^r (A) P:,(2A + tp) K^'iX + ^p)} 

A[;^*+(A)] = ^Ta,{P--,KL{X + ^p)Rl,{-2X-3^p)KL{X)]. (A.ll) 

One obtains similar relations when the spaces a and b are fused. To compute the quantum determi- 
nants explicitly we use the following identities which can be easily proved with the help of unitarity 
()1.5j) and the crossing relation p.9j) 

P--, R-aM) R,M + W) = C(A + zp) P--, , (A.12) 

P;^ RaM) R-UX + ip) = C(A + ip) P;^, m = l,...,L* (A.13) 
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where L* = L in the SP case and L* = 2L in the SNP case. One then writes when fusing the spaces 
a and b 

6[T{\)] = C(A + 2p)^*/2r(A + ^p)^*/2 , S[f*{\)] = C(A + ^p)'^*/\{\ + ipf*" (A.14) 

whilst, when we fuse the spaces a and 6, 

5[T{\)] = C(A + 2p)^*/'C(A + tpf'^ , 5[r*(A)] = C(A + zp)^*/'C(A + tpf*'^ ■ (A.15) 

Furthermore, the role of C, and C is interchanged in the latter equation depending whether the space 
Vm belongs to the fundamental representation or to its conjugate. This statement is important if one 
aims at constructing the alternating spin chain. Finally for the special case K~ = 1 and K^ = 1 

A[K*-{X)] = q*{2X + ip), A[K*+{X)] = q*{-2X~3ip) , (A.16) 

where 

g*(A) = q{X) for SP, g*(A) = g(A) for SNP 
g(A) = A-ip, q{X) = X + i. (A.17) 

B Generalised fusion 

We describe a generalised fusion procedure for sl{J\f) open spin chains [42]. The procedure we 
use follows the lines of the construction of the Sklyanin determinant for twisted Yangians [37] and 
reflection algebras [40,52]. The crucial observation here is that for the general case an one dimensional 
projector can be also obtained by repeating the fusion procedure M times, this is because J^f^-^ = 
1 © ... . The procedure described in the previous section is basically consequence of the fact that 
Ar®A/'= 1 © (AT^-l). 

Let us now introduce the following necessary objects for the generalised fusion procedure for open 
spin chains (see also equations (2.13), (2.14) in [53]), 

T<a> = T<^ai...a^> = ^^^(Ai) . . .Ta^{Xj^), T*„> = T*^{Xi) . . .T„*^(Aa/-) (B.l) 

where A; = A + z(/ — 1), I = 1, . . . ,Af and R* defined in flA.l|) . For two sets {p/}«=i,...,A/' and {pJ}/=i,...,Ar 
we also define 

k=l,...,Af-l 

In particular, 7^<a>({Ai}, {A;}) is proportional to the antisymmetriser A, i.e. the projector onto a 
one-dimensional space {A^ = A) . We also use A'^ = I — A and 

T^*<a>{M) = n R:,aj,i-f^k-PM-2ip)...Rl,^J-pk-Pk+i-2tp) (B.3) 

fe=l,...,A/'-l 
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By multiplying the four equations 



-^ <a> ^ vA i <a.> + vA i <a> , 

f ^,> = ^ f *„> + ^+ f *„> (B.4) 

and keeping in mind that 

A 7^*<„> A+ = A 7^*<+> ^+ = ^ T<„> ^+ = ^ f^,> ^+ = (B.5) 

we get 

Tr<„>(7^^+> T<„> 7^^,> f *„>) = Tr<„>(^ 7^*<+> ^ T<„> ^ 7^^,> A f *„>) 

+ Tr<„>(^+ 7^^+> ^+ T<„> ^+ 7^*<„> X f *„>). (B.6) 

Applying equation 

T,(AJ R:,{Xa + A,) f;(A,) = f;(\) i?:,(A, + \) T„(AJ (B.7) 

recursively we can show that 

T<a> 7^*<„> f^„> = T<*,> (B.8) 

where 

^<a>= n P«:(^'^) n ^;a.(A. + A.)) . (B.9) 

A:=l,...,Ar V l=k+l,...,M J 

However, as discussed in [53] the trace of the above quantity decouples to a product of M trans- 
fer matrices, and therefore the LHS of ()B.6j) simply becomes 11/=! ^(Ai)- Taking into account the 
property, 

A 0<a> A = Tr<«>(^ 0<a>) A, (B.IO) 



we can write the first term of the RHS of ()B.6jl as product of quantum determinants, which are 
simply c numbers, i.e. 

Tr<„>M Ka> ^ T^a> A 7^^„> A t*„>) = A{K*+(A)} 6{T{X)} A{K*-{X)} 6{f*iX)} (B.ll) 

where 

A{ir*+(A)} = Tr<a>M7^*<+>}, 6{TiX)} = Tr^a>{AT^a>}, 

A{K*-{X)} = Tr<„>M7^*<„>}, (5{t*(A)} = Tr<„>M f^^J. (B.12) 



Finally the second term of the RHS of ()B.6|) is simply the fused transfer matrix t(A). Therefore, 
equation ()B.6|) can be rewritten as 

iW = n^(^') - ^{^^*^(A)} HT{X)} A{/r-(A)} 5{f*{X)}. (B.13) 



1=1 



31 



References 

[1] R.J. Baxter, Partition function of the eight-vertex lattice model, Ann. Phys. 70 (1972) 193; 
J. Stat. Phys. 8 (1973) 25; 
Exactly solved models in statistical mechanics (Academic Press, 1982). 

[2] L.D. Faddeev, N.Yu. Reshetikhin and L.A. Takhtajan, Quantization of Lie groups and Lie 
algebras, Leningrad Math. J. 1 (1990) 193. 

[3] I.V. Cherednik, Factorizing particles on a half line and root systems, Theor. Math. Phys. 61 
(1984) 977. 

[4] E.K. Sklyanin, Boundary conditions for integrable quantum systems, J. Phys. A21 (1988) 2375. 

[5] P.P. Kuhsh and E.K. Sklyanin, Algebraic structure related to the reflexion equation, J. Phys. 
A25 (1992) 5963. 

[6] P.P. Kuhsh and R. Sasaki, Covariance properties of reflection equation algebras, Prog. Theor. 
Phys. 89 (1993) 741 and hep-th/9212007. 

[7] J. Donin, P. P. Kuhsh and A. I. Mudrov, On universal solution to reflection equation, Lett. 
Math. Phys. 63 (2003) 179 and math. QA/0210242. 

[8] H.J. de Vega and A. Gonzalez-Ruiz, Boundary K-matrices for the XYZ, XXZ and XXX spin 
chains, J. Phys. A27 (1994) 6129 and hep-th/9306089. 

[9] S. Ghoshal and A.B. Zamolodchikov, Boundary S matrix and boundary state in two-dimensional 
integrable quantum field theory. Int. Journ. Mod. Phys. A9 (1994) 3841. 

[10] J. Abad and M. Rios, Nondiagonal solutions to reflection equations in SU{N) spin chains, Phys. 
Lett. B352 (1995) 92 and hep-th/9502129. 

[11] A. Doikou and R.I. Nepomechie, Bulk and Boundary S Matrices for the SU(N) Chain, Nucl. 
Phys. B521 (1998) 547 and hep-th/9803118. 

[12] P. Bowcock, E. Corrigan, P.E. Dorey and R.H. Rietdijk, Classically integrable boundary condi- 
tions for afEne Toda field theories, Nucl. Phys. B445 (1995) 469 and hep-th/9501098; 
P. Bowcock, E. Corrigan and R.H. Rietdijk, Background field boundary conditions for affine 
Toda field theories, Nucl. Phys. B465 (1996) 350 and hep-th/9510071. 

[13] G.M. Gandenberger, New non- diagonal solutions to the An boundary Yang-Baxter equation, 
hep-th/9911178. 

[14] A. Doikou, Quantum spin chain with "soliton nonpreserving" boundary conditions, J. Phys. 
ASS (2000) 8797 and hep-th/0006197. 



32 



[15] H. J. de Vega and A. Gonzalez- Ruiz, Exact solution of the SUq{n) invariant quantum spin 
chains, Nucl. Phys. B417 (1994) 553 and hep-th/9309022. 

[16] A. Doikou, Fusion and analytical Bethe Ansatz for the A„_i open spin chain, J. Phys. A33 
(2000) 4755 and hep-th/0006081. 

[17] H.J. de Vega and F. Woynarovich, New Integrable Quantum Chains combining different kind of 
spins, J. Phys. A25 (1992) 4499. 

[18] H. Saleur and B. Wehefritz-Kaufmann, Integrable quantum field theories with 0SP{m,\2n) sym- 
metries, Nucl. Phys. B628 (2002) 407 and hep-th/0 112095. 

[19] D. Arnaudon, J. Avan, N. Crampe, A. Doikou, L. Frappat and E. Ragoucy, Classification of 
reflection matrices related to (super) Yangians and application to open spin chain models, Nucl. 
Phys. B668 (2003) 469 and math.QA/0304150. 

[20] D. Arnaudon, J. Avan, N. Crampe, A. Doikou, L. Frappat and E. Ragoucy, Bethe Ansatz 
equations and exact S matrices for the osp{M\2n) open super spin chain, Nucl. Phys. B687 
(2004) 257 and math-ph/0310042. 

[21] A. Gonzalez- Ruiz, Integrable open-boundary conditions for the supersymmetric t-J model. The 
quantum group invariant case, Nucl. Phys. B424 (1994) 468 and hep-th/9401118. 

[22] F.H.L. Essler, The Supersymmetric t — J Model with a Boundary, J. Phys A29 (1996) 6183 
and cond-mat/9605180. 

F.H.L. Essler, V.E. Korepin, Higher conservation laws and algebraic Bethe Ansatze for the 
supersymmetric t — J model, Phys. Rev. B46 (1992) 9147. 

F.H.L. Essler, V.E. Korepin, A New Solution of the Supersymmetric t — J Model by Means of 
the Quantum Inverse Scattering Method, and hep-th/9207007. 

[23] B.Y. Hou and R.H. Yue, General solution of reflection equation for eight vertex model, Phys. 
Lett. A183 (1993) 169. 

B.Y. Hou, H. Fan and R.H. Yue, Exact diagonalization of the quantum supersymmetric 
SUg{n\m) model, Nucl. Phys. B462 (1996) 167. 

T. Deguchi and R.H. Yue, Magnetic susceptibility and low-temperature specific heat of inte- 
grable 1-D Hubbard model under open boundary conditions, J. Phys. A30 (1997) 8129. 
T. Deguchi, K. Kusakabe and R.H. Yue, A gapless mode induced by the boundary state in the 
half-filled Hubbard model, J. Phys. A31 (1998) 7315. 

B.Y. Hou, R.H. Yue and X.Q. Xi, Exact solution of 1-D hubbard model with impurities under 
the open boundary condition, Phys. Lett. A257 (1999) 189. 

[24] A. B. Zamolodchikov and Al. B. Zamolodchikov, Factorized S-matrices in two dimensions as 
the exact solutions of certain relativistic quantum field theory models. Annals Phys. 120 (1979) 
253. 



33 



[25] M. Gattobigio, A. Liguori and M. Mintchev, Quantization of the nonlinear Schrodinger equation 
on the halfhne, Phys. Lett. B428 (1998) 143 and hep-th/9801094; 

The nonhnear Schrodinger equation on the half line, J. Math. Phys. 40 (1999) 2949 and 
hep-th/9811188. 

[26] N.J. MacKay and B.J. Short, Boundary scattering, symmetric spaces and the principal chiral 
model on the half-line, Commun. Math. Phys. 233 (2003) 313 and hep-th/0104212. 

[27] L. Mezincescu and R.I. Nepomechie, Analytical Bethe Ansatz for quantum algebra invariant 
spin chains, NucL Phys. B372 (1992) 597 and hep-th/9110050. 

[28] V.G. Drinfel'd, Hopf algebras and the quantum Yang-Baxter equation, Soviet. Math. DokL 32 
(1985) 254; 
A new realization of Yangians and quantized afhne algebras, Soviet. Math. DokL 36 (1988) 212. 

[29] H. Saleur, The continuum limit of SL{N\K) integrable super spin chains, NucL Phys. B578 
(2000) 552 and solv-int/9905007. 

[30] M. J. Martins, Integrable mixed vertex models from braid-monoid algebra. Statistical Physics on 
the Eve of the 21st century, Ed. M.T. Batchelor and L.T. Wille, Series on Advances in Statistical 
Mechanics, vol 14, World Scientific, Feb 1999, and solv-int/9903006. 

[31] J.B. McGuire, Study of exactly soluble one-dimensional N-body problems, J. Math. Phys. 5 
(1964) 622. 

[32] C.N. Yang, Some exact results for the many-body problem in one dimension with repulsive delta- 
function interaction, Rev. Lett. 19 (1967) 1312. 

[33] V.E. Korepin, New effects in the massive Thirring model: repulsive case, Comm. Math. Phys. 
76 (1980) 165; 

V.E. Korepin, G. Izergin and N.M. Bogoliubov, Quantum inverse scattering method, correlation 
functions and algebraic Bethe Ansatz (Cambridge University Press, 1993). 

[34] L.D. Faddeev and L.A. Takhtajan, Spectrum and scattering of excitations in the one-dimensional 
isotropic Heisenberg model, J. Sov. Math. 24 (1984) 241; 

L.D. Faddeev and L.A. Takhtajan, What is the spin of a spin wave? Phys. Lett. A85 (1981) 
375. 

[35] M. Mintchev, E. Ragoucy and P. Sorba, Spontaneous symmetry breaking in the gl(N)-NLS 
hierarchy on the half line, J. Phys. A34 (2001) 8345 and hep-th/0104079. 

[36] G.I. Olshanski, Extension of the algebra tl{Q) for infinite dimensional classical Lie algebras Q, 
and the Yangians U{gl{m)), Soviet. Math. DokL 36 (1988) 569; 

Representations of infinite dimensional classical groups, limits enveloping algebras and Yangians, 
in "Topics in Representation theory", A. A. Kirillov ed.. Advances in Soviet. Math. 2 (1998) 1. 

34 



[37] A. Molev, M. Nazarov and G. Olshanski, Yangians and classical Lie algebras, Russian Math. 
Survey 51 (1996) 205 and hep-th/9409025. 

[38] D. Arnaudon, R. Poghossian, A. Sedrakyan and P. Sorba, Integmhle chain model with additional 
staggered model parameter, Nucl. Phys. B588 (2000) 638 and hep-th/0002123 
D. Arnaudon, A. Sedrakyan and T. Sedrakyan, Multi-leg integrable ladder models, Nucl. Phys. 
B676 (2004) 615 and h ep-th/0210087, 

[39] H.J. de Vega and A. Gonzalez- Ruiz, Boundary K-matrices for the six vertex and the n{2n — 1) 
An-i vertex models, J.Phys. A26 (1993) L519 and hep-th/9211114. 

[40] A. I. Molev and E. Ragoucy, Representations of reflection algebras, Rev. Math. Phys. 14 (2002) 
317 and math.QA/0107213. 

[41] D. Arnaudon, J. Avan, N. Crampe, L. Frappat and E. Ragoucy, R-matrix presentation for 
super-Yangians Y {osp{m\2n)) , J. Math. Phys. 44 (2003) 302 and math ■QA/0111325. 

[42] N. Crampe, Approches algebriques dans les systemes integrables. University of 
Savoie PhD Thesis (June 2004), in French, Preprint LAPTH-These-1056/04, 
http : //wwwlapp . in2p3 . f r/preplapp/psth/LAPTH1056 . ps . gz. 

[43] W. Galleas and M.J. Martins, Solution of su{N) vertex model with non-diagonal open bound- 
aries, nlin. SI/0407027. 

[44] V.G. Kac, Lie super algebras. Adv. Math. 26 (1977) 8; A sketch of Lie superalgebra theory, 
Commun. Math. Phys. 53 (1977) 31. 

[45] D.A. Leites, M.V. Saveliev, V.V. Serganova, Embeddings of Lie superalgebra osp{l\2) and the 
associated nonlinear supersymmetric equations, in Group theoretical methods in physics, M.A. 
Markov, V.I. Man'ko and V.V. Dodonov (eds.), p. 255, VUN Science Press, Utrecht, The Nether- 
lands (1986). 

[46] V.K. Dobrev and V.B. Petkova, On the group-theoretical approach to extended conformal su- 
persymmetry : function space realizations and invariant differential operators, Fortschr. d. Phys. 
35, 537-572 (1987). 

[47] L. Frappat, A. Sciarrino, P. Sorba, Dictionary on Lie algebras and super algebras, Academic 
Press (London), 2000. 

[48] C. Briot and E. Ragoucy, Twisted superYangians and their representations, J. Math. Phys. 44 
(2003) 1252 and math. QA/01 11308. 

[49] A. Foerster and M. Karowski, The supersymmetric t — J model with quantum group invariance, 
Nucl. Phys. B408 (1993) 512; 

A. Foerster, The supersymmetric t — J model and its quantum group deformation, PhD Thesis, 
Berlin University, April 1993. 

35 



[50] N. Yu. Reshetikhin, The spectrum of the transfer matrices connected with Kac-Moody algebras, 
Lett. Math. Phys. 14 (1987) 235. 

[51] L. Mezincescu and R.I. Nepomechie, Fusion procedure for open chains, J. Phys. A25 (1992) 
2533. 

[52] A. I. Molev, Yangians and their apphcations, in "Handbook of Algebra", VoL 3, (M. Hazewinkel, 
Ed.), Elsevier, 2003, 907, math-QA/0211288. 

[53] J. Avan and A. Doikou, Commuting quantum traces: the case of reflection algebras, J. Phys. 
A37 (2004) 1603 and math.QA/0305424. 



36 



